ON PARTICLE TRAJECTORIES IN LINEAR DEEP-WATER WAVES 



ANCA-VOICHITA MATIOC 

Abstract. We determine the phase portrait of a Hamiltonian system of equations describing 
the motion of the particles in linear deep-water waves. The particles experience in each period 
a forward drift which decreases with greater depth. 



1. Introduction 

The motion of water particles under the waves which advance across the water is a classical 
problem in this field. Watching the sea it is oft possible to trace a wave as it propagates on 
the water's surface, but what one observes traveling across the sea is not the water but a wave 
pattern. Although the wave travels from one place to another, the substance through which 
it travels moves very little. As the wave advances across the water and can be followed for a 
long way, a typical water particle moves slightly up and down, forward and backward as the wave 
passes it. If an object hovered in the water, like a water particle, its motion would be synchronized 
with that of a floating object lying on the water's surface, with its orbit diminishing with the 
distance from the surface. As waves generated by wind in an area move towards a region where 
the wind has ceased, we observe swell-long crested two-dimensional waves approaching a smooth 
sinusoidal shape and moving over long distances. Deep-water waves are modelled mathematically 
as periodic two-dimensional waves in water of infinite depth. The motion of the water particle 
in the fluid below swell is of great interest. The classical description of these particle paths is 
obtained within the framework of linear water wave theory [1, 15, 16, 23, 24, 25, 29, 30]: all 
water particles trace a circular orbit, the diameter of which decreases with depth so that the 
orbital motion practically ceases at depth equal to one-half the wavelength. These features have 
important practical consequences. For example, a submarine at a depth below half a wavelength 
would hardly notice the motion of the surface wave, for this reason submarines dive during storms 
in the open sea. 

The only known explicit solution with a non-flat free surface of the governing equations for 
gravity water waves is Gerstner's wave [17]: a deep-water wave solution for which all particle paths 
are circles of diameters decreasing with the distance from the free surface (see the discussion in 
[2, 3]). Due to the mathematical intractability of the governing equations for water waves, for 
irrotational water waves (Gerstner's wave has a peculiar nonvanishing vorticity) the classical 
approach [15, 20, 23, 24, 28, 29, 30] relies on analyzing the particle motion after linearization 
of the governing equations. However, even within the linear water wave theory, the ordinary 
differential equations system describing the motion of the particles is nevertheless nonlinear and 
explicit solutions of this system are not available, but qualitative features of the underlying flow 
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are known [6, 9, 10, 11, 21, 26, 27, 33] and could possibly lead to a confirmation of the features 
we prove here, within the confines of linear water waves (see [4]). 

In this paper we show that for linear deep-water wave no particle trajectory is actually closed, 
unless the free surface is flat. Each trajectory involves over a period a backward/forward move- 
ment of the particle, and the path is an elliptical arc with a forward drift. Support for this 
conclusion is given by the analysis of the average flow of energy within linear water wave theory 
(see [23]): due to the passage of a periodic surface wave, the water particles in the fluid experience 
on average a net displacement in the direction in which the waves are propagating, termed Stokes 
drift [31] (see the discussion in [34]). This result was first obtained in [13] for gravity waves and 
[19] for capillary and capillary-gravity water waves. This results have been improved in [4, 9, 18] 
in the case of irrotational waves by describing the exact geometry of the actual particle paths. For 
gravity deep-water waves we reconfirm herein the results of [13], but our analysis is more precise 
and uses only elementary analysis methods. The work is structured as follows: we recall first 
the governing equations for water waves and give the exact solutions of the linearised problem. 
Our major contribution, in Section 3, is a precise analysis of the phase portrait of a Hamiltonian 
system associated to the particle motion. This is the key point in determining the trajectories of 
the water particles. 

2. Preliminaries 

In this section we recall the governing equations for the propagation of two-dimensional gravity 
deep-water waves and we present their linearization (for a more detailed discussion we refer to 
[23]). 

2.1. The governing equations. We consider a two-dimensional inviscid incompressible fluid in a 
constant gravitational field. To describe the waves propagation we consider a cross section of the 
flow that is perpendicular to the crest line with Cartesian coordinates (x,y), the y-axis pointing 
vertically upwards and the x-axis being the direction of wave propagation (see Figure 1). Let 
(u(t,x,y),v(t,x,y)) be the velocity field and the equation of the free surface with mean water 
level zero is y = rj(t,x) with J„rj(t,x) = 0. 

Constant density (homogeneity) is a physically reasonable assumption for gravity waves [23], 
and it implies the equation of mass conservation 

u x + u y = 0. (2.1) 

Under the assumption of inviscid flow the equation of motion is Euler's equation 

U t + UU X + VUy = -P x 

(2.2) 

v t + uv x + Wy = -Py - g, 

where P(t, x, y) denotes the pressure and g is the gravitational constant of acceleration. The free 
surface decouples the motion of the water from that of the air, a fact that is expressed in [23] by 
the dynamic boundary condition 

P = P on y = V (t,x), (2.3) 

if we neglect surface tension, where Pq is the (constant) atmospheric pressure. Since the same 
particles always form the free surface, we also have the kinematic boundary condition 

v = T] t + ur] x on y = r](t, x). (2.4) 
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FIGURE 1. A deep-water wave with wavelength A 



The boundary condition at the bottom 

(u,v) —> (0,0) as 1/ — > -co uniformly for i£K, (2-5) 

expresses the fact that at great depths there is practically no motion. Summarising, the governing 
equations for water waves are encompassed by the nonlinear free-boundary problem 
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(2.6) 



An important category of flows are those of zero vorticity (irrotational flows), characterized by 
the additional assumption 

u y = v x . (2.7) 
The vorticity of a flow, u = u y — v x , measures the local spin or rotation of a fluid element so that 
in irrotational flows the local whirl is completely absent. Relation (2.7) ensures the existence of 
a velocity potential cj)(t, x, y) defined up to a constant via 

(j) x =U, (j)y= V. 

In view of (2.1), is a harmonic function, i.e. (d x + dy)(j) = so that the methods of complex 
analysis become available for the study of irrotational flows. 

Concerning the physical relevance of irrotational water flows, field evidence indicates that for 
waves entering a region of still water the assumption of irrotational flow is realistic [25]. Moreover, 
as a consequence of Kelvin's circulation theorem [1], a water flow that is irrotational initially 
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has to be irrotational at all later times. It is thus reasonable to consider that water motions 
starting from rest will remain irrotational at later times. Nonzero vorticity is the hallmark of the 
interaction between waves and non- uniform currents cf. [10] and [12, 22]. The most important 
example of rotational flows are tidal currents cf. the discussion in [5, 32], but the assumption of 
constant vorticity can not be accomodated within the framework of deep-water flows. However, 
we consider herein irrotational gravity waves. 

2.2. Linear water waves. The problem (2.6) is nondimensionalised using a typical wavelength A 
and a typical amplitude of the wave e. We define the set of nondimensional variables 



A 

t*-> -=t, 
V9 
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x h-» Ax, y i— > y, 
combined with the scaling 

pt-^ep, (u, v) i— > e(u, v). 
Setting the constant water density p = 1, the pressure in the new nondimensional variables is 

P = Pq- gy + gp, 

with the nondimensional pressure variable p measuring the deviation from the hydrostatic pressure 
distribution. We obtain the following boundary value problem in nondimensional variables 
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The linearized problem is now obtained by letting e — > in (2.8). In the limit we arrive at 
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(2.9) 



v = Vt 

(u,v) — > (0,0) as y — > — oo. 

It is worth noticing that the regularity results on the streamlines for water waves obtained in the 
recent papers [27] and [8] ensure that good approximate properties hold. 

Looking for solutions of (2.9) representing waves traveling at speed c > 0, we impose that all 
functions u,v,p and ij have an (x, independence in the form of x — ct. Furthermore, one seeks 
even periodic functions of period one, the evenness expressing the symmetry of the surface wave: 
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it is known, cf. [6] and [7, 14] that a gravity surface traveling wave whose profile is monotone 
between the crests and troughs has to be symmetric about the crest. Thus, we are led to the 
fundamental Fourier mode Ansatz 

r](x,t) = cos[27r(x — ct)]. 
For this specific rj the problem (2.9) has the solution 

f f]{t,x) = cos[27r(x — ct)] 
u(t,x,y) = Ap. cos[27r(x — ct)], 



v(t,x,y) 
P{t,x,y) 



f(y) sin[2vr(x - ct)], 



A 



cos[27r(x — ct)], 




whereby the constant c and the function / are given by 

/27T 

7^, f(y) = 2vrcexp ( —y 

Marking the variables of this nondimensional solution by tildes, we recover the physical vari- 
ables by the backward transformation 

A ~ 

x = \x, y = y, t = —t, 

V9 



n ~ Va- 
Defining the wavenumber k and the frequency u by 



7] = er], p = ep. 



k 



2vr 
T' 
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gk, 



where 2ir(x — ct) = kx — ut, we finally obtain that the tupel 

rj(t,x) = e cos(kx — cut), 

u(t,x,y) = euj exp(ky) cos(kx — oot) , 

v{t,x,y) = eco exp(ky) sm(kx — cot) , 

P(t,x,y) = P - gy + egexp(ky)cos(kx - ujt) 



(2.10) 



is a solution of the linearised deep-water wave problem. Note that in the physical variables the 
speed of the wave in (2.10) is 
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3. Particle trajectories 

In this last section we study the trajectories of the water particles in the linear wave (2.10). 
First of all, we note that a fluid particle (x(t),y(t)) must satisfy the following equations 

dx dy 

so that, in view of (2.10), the motion of the particle is described by the system 

M exp(/cy) cos(kx — ut), 

(3.1) 

M exp(/cy) sin(/cx — tot), 

with initial data (x(0),y(0)) = (xo,yo)- We use herein the shorthand 

M := eu. (3.2) 

The right-hand side of the differential system (3.1) is smooth so that the existence of a unique 
local smooth solution is ensured the Picard-Lindelof theorem. Without actually solving (3.1) we 
would like to display the principal features of the solution. To study the exact solutions to (3.1) 
it is convenient to re-write the system in a moving, and we define therefore the variables 

X = kx-ut, Y = ky. (3.3) 

With this notation (3.1) is equivalent to the following problem 

dX 

— = kM exp(Y) cos(X) - u, 

dY (3 - 4) 

— = A;Mexp(y)sin(X), 

with (X(0),Y(0)) = k(xo,yo). In order to determine the phase portrait corresponding to (3.4), 
we take into consideration the property of (3.4) of being a Hamiltonian system. Indeed, setting 

H(X, Y) = kM exp(Y) cos(X) - uY, (X, Y) £ M 2 , (3.5) 

we re-express problem (3.4) as follows 

^ Q TT 

f (3.6) 
f - 

Since H is constant on solutions of (3.6), we are left to determine the level curves of the function 
H. In [13] the authors studied the phase portrait of (3.6) approximately, by using Morse's lemma. 
A precise analysis allows us to determine herein the phase portrait of the Hamiltonian system 
of equations describing the motion of the particles in linear deep-water waves by using only 
elementary methods. 

Note that H is 2tt— periodic in X and even with respect to the same variable. Therefore, 
the level curves of H are determined by the level sets of the restriction H : 0, — > R, where 
:= [0,7r] x R. Furthermore, it is not difficult to see that the equation Vif = has a unique 
solution 

P:=(X„y,):= (0>(^ 
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FIGURE 2. Phase portrait in the moving frame (a\ < a* < 0.2) 

meaning that P is the unique stationary solution of (3.6) within 0. Moreover, P is located on 
the level curve -ff _1 ({a*}) with a* given by 



u 1 



\kMJJ 



We make now a crucial observation. Let a S M be given with H(X, Y) = a for some (X, Y) £ Q. 
In view of fcMexp(Y~) cos(X) — uY = a, it must hold that 

X = f a {Y) :=arccos ( fc ^[ y ) 

with / a : -D(/ Q ) C R — >■ [—1,1]. The definition domain D(f a ) of f a consists of all the points 
Y 6 R which satisfy the inequalities 

a + uY 



1 < 9a{Y) 



< 1. 



(3.7) 



£;Mexp(Y) 

Hence if" 1 ({a}) flfi = graph for all a 6 R, and we are left to determine the definition 
domain of the function f a , that is to solve the inequality (3.7). The first important result of this 
paper is the following theorem which gives a precise description of the phase curves of (3.6): 

Theorem 3.1. Given a € R, there exists a function f a : D(f a ) Cl->l such that if _1 ({a}) = 
graph(/ a ) for all a € R. More precisely, there exists a unique Y n (a) > such that 
(a) If a< a*, then D(f a ) = [Y n (a),oo), 

f a {Y-K{a)) = tt, and 



Jim f a {Y) = tt/2. 



Moreover, the minimum of f a belongs to the interval (0, 7r/2). 
(6) If a > a*, there exist positive real numbers Yi(a), i = 1,2, l^-(a) < ^i(ct) < ^2(0) 
suc/t i/iaf D(/ Q ) = [Y^aOjYita)] U [12(a), 00). 
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Moreover, f a : [Y n (a),Yi(a)] — > [0,tt] is bijective and decreasing, whereas f a : [Y2(a),oo) — > 
[0, vr/2) is bijective and increasing. 
Additionally, we have 

(1) Y„- : M — ?■ R is bijective and decreasing; 

(2) Y\ : (a*,oo) — > (— oo,Y.) is bijective and decreasing; 

(3) Y2 : (a*, 00) —7- (Y,,oo) is bijective and increasing; 

(4) {Yi{a) - Y n (a)) \0asa/oo. 

The functions mentioned above are real analytic in the interior of their definition domain (ex- 
cepting f at which is real analytic only on (Y^a*), Y.) U (Y*,oo)). 

Remark 3.2. Having proved Theorem 3.1, the phase portrait is complete (see Figure 2). Note 
that by (4), we know that the phase curves parametrised by f a : [Yr(a), Yi(a)] — > [0,7r] flatten 
out as q 00. The motion along the phase curves is determined by the fact that 

-d x H{X,Y) = kMexp(Y) sm(X) > 

provided X G (0, it). The phase portrait also discloses that the stationary solution (X*, Y*) is a 
saddle node for (3.5). Solutions starting in a point (Xq,Yq) located on a separatrix (one of the 
level curves containing P) are attracted by (X*,Y*) provided Yo < Y*. On the other hand, if 
Yq > Y* they tend to infinity. 

Proof of Theorem 3.1. Let a G K. We observe first that the equation g a (Y) = — 1 is always 
solvable in IR. Indeed, it is equivalent to the following relation 

a = h{Y) := -Mkexp(Y) - ojY. 

Clearly, the function h : M — > M is strictly decreasing and, moreover, 

lim h(y) = 00 and lim h(y) = —00. (3-8) 
Y— >— 00 Y— >oo 

Whence, there exists a unique real number Yr(a) such that 5 a (Y,-(a)) = — 1. Consequently, all 
phase curves of (3.6) intersect the line [X = ir] exactly once, that is at 

ff -1 ({a})nLY = 7r] = Y 7r (a) for all a G R. 

A simple calculation shows that 

' (Y ( )) - — !— ( - a + u}Y n (a) \ _ w + KM exp(Y 7r (q)) 

9 a { *(<*)) - kM ^ exp(y ^ (a)) exp(Y w (a)) ) ~ kM exp(Y^(a)) ' 

and we infer from the implicit function theorem that Y n is real analytic and, additionally 

+ ^(Y 7r (a))Y» = 0. (3.9) 



kM exp(Y 7r (a)) 

This implies that Y v is strictly decreasing. Due to (3.8), we conclude that Y n is bijective. Differ- 
entiating g a at Y G M, we observe that 

^-a-^Y f <0, if Y>(^-a)/ W , 
5al J fcMexp(Y) 1 > , if y<(w-a)/w. 1 J 
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Let us now compare the value of g a at the threshold value (u — a) /id, where the monotonicity 
changes, to 1, the larges value g a can take in order to remain in the definition domain of arccos. 
We recover a*, since 

g a ((id — a)/uj) < 1 if and only if a < a*. 

Analogously, we have that g a ((id — a) /id) > 1 if and only if a > a*. Since g a ((id — a) /id) is the 
(positive) maximum of g a and limy_ 5 ._ 00 g a (Y) = — oo, we conclude that 

Y^(a) < (u — a)/u 

for all a£l. Particularly, if a < a* then g a (Y) < 1 for all Y > Y n (a). Moreover, taking into 
account that limy^oo g a (Y) = 0, we conclude that D(f a ) = [Y n (a), oo) for all a < a*. 

On the other hand, if a > a*, then g a {{oj — ot)/oS) > 1, and there exist precisely two points 
Y\{q) < (d — a)/u < 12(a) with the property that g a (Yi(a)) = 1, i = 1,2. By (3.10) we have 
that g' a {Y\{a)) > 0, g' a (Y2(a)) < 0, and we infer from the implicit function theorem that Y\ is 
strictly decreasing, while Y2 is strictly increasing. We have use the fact that both Y\ and Y2 
verify, due to g a (Yi(a)) = 1, i = 1, 2, relation (3.9). With this, the assertions (1) and (2) are 
immediate. 

We finish the proof, by noticing that for a sufficiently large 

2 = a (Yi(a)) - g a (Y v (a)) = ^ a (y c (a))(Yi(a) - Y w (a)). 
for some Y c (a) G (l^-(a), Yi(a)). Consequently, we have 

w — a — wli (a) ,, , . . . . . 
2£ lM^W (Fl(Q) " F - (Q)) 

and, for a — )■ 00, we obtain in view of Theorem 3.1 (3) that Y\{a) — Y n (a) — > a -^oo 0. We show 
now that the convergence is monotonia Indeed, from (3.9) and (3.10) we see that 

K( a ) = 1 w n and Y (( a ) 



uj — a — ujY n (a) oj — a — ajY\(a) 

Consequently, 

w(y ff (a)-Yi(a)) 



Y{(a)-Y;(a) 



(id — a — (jjY K (a)) (oj — a — wYi(a)) 
2cj 



< 0. 



(w — a — idY n (a)) (id — a — idY\(a)) 
This completes the proof. □ 

Knowing the phase curves (X(t),Y(t)), the particle trajectories (x(t),y(t)) in the linear wave 
(2.10) are given, via (3.3), by 

X(t)= k + k f > V{t)= k- (3 ' U) 

For physical reasons realistic particle paths are bounded. Thus, the critical point P must at least 
lie above the typical height of the wave. Since the mean surface level is zero and the amplitude is 
e, a crest has y— coordinate e. Since the y— coordinate of P is ln(l/A:£:)/A:, we must thus impose 
that ek < ln(l/efc), that is 

ekexp(ek) < 1. (3.12) 
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FIGURE 3. Particle trajectory above the flat bed 

This is the only condition which has to be imposed for (2.10) to be a realistic model of water 
waves of infinite depth, wave number k and small amplitude e. This is an improvement compared 
to [13], because we do not have to to make sure that the separatrix intersects the the line X = n 
only once; this is in our case a conclusion of studying the phase portrait of (3.6). 

We describe now the trajectory of the water particle (x, y) in the fluid. We may consider the 
orbit of a particle (xo>yo) which starts at the lowest point (corresponding to the line X = ir in 
the phase-portrait). This trajectory corresponds to the solution (X, Y) of (3.6) which starts at 
(Xq,Yq) = (kxo,kyo)- Furthermore, we denote by 9 = 6(Yq) the time needed for (X, Y) to reach 
the line X = —ir. Indeed, from thre first equation of (3.6) we have for Y located beneath Y ir (a^), 
that is under the separatrix, that 

r ds 

9{Y) = 2 J ^-feMexp(y)cos( S )' (3 ' 13) 

hence 9 depends only of Y . Let us observe the evolution of the particle (xo, Vo) as t G [0, 0\. Since 
X = kx — out, Y = ky, we find see from (3.6) that 

^<0, ^<0 forX(t)€(-7r,-7r/2), 
^>0, ^<0 forX(i)e(-7r/2,0), 
^|>0, ^>0 forX(i)e(0,7r/2), 

^<0, ^>0 forX(t)€(7r/2,7r). 

Also, from the phase portrait of (3.6), we see that the hight y of the particle reaches its initial 
value for the first time when t = 9. What we still have to establish is the sign of x(9) — x(0). 
Since 

x{9) - X ( ) = 9 -^L > o, 

k 

we find from (3.13) that 

_ , (0) _ l f t Mex P (F)co 3 (») ^ _ 2 r cos(s) & 4 ^ s ^(x) 



k J to — kM exp(y) cos(s) k J a — cos(s) k J a 2 — sin 2 (x)' 
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with a := U)/{kM exp(Y)). In view of X' < 0, we conclude that the water particle experiences a 
forward drift equal to x(6) — x(0) > 0. 

Clearly, by the definition of a we see that the drift decreases with depth and, in the limit 
Y — > — oo, the forward drift vanishes x{9) — x(0) —> 0. 

Whence, we have shown that: 

Theorem 3.3. The particle paths in the linear water wave (2.10) are not closed. There exists a 
forward drift over a period, and this drift decreases with depth. 
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